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1. Introduction 

Recently, a new approach to the string representation of the confining phase of gauge theories 
was proposed [1]. This is the so-called Universal Confining String Theory (UCST), which is the 
theory of an Abelian antisymmetric tensor field with a nonlinear action, interacting with the string. 
This field effectively substitutes an infinite number of monopoles in the 3D compact QED. It was 
argued in Ref. [1] that the summation over the branches of the UCST action should correspond 
to the summation over the string world-sheets. It was also proved in Ref. [1] that the UCST 
partition function, which is nothing else but the Wilson average in the 3D compact QED, satisfies 
loop equations [2] modulo contact terms. However, the statement made in Ref. [1] concerning 
the universality of the wave operator standing on the L.H.S. of the usual loop equations obtained 
from the Yang-Mills theory, which were derived and investigated in Ref. [2], was only a conjecture, 
and the relation of these equations to the equation of motion of the tensor field in the UCST was 
absolutely unclear. In order to clarify this conjecture, in Ref. [3] the loop equation for the 4D 
UCST partition function was derived and investigated. In particular, it was demonstrated in [3] 
that the wave operator of the obtained loop equation is quite different from the usual one. Also 
the corresponding contact term was calculated explicitly. 

The low-energy limit of the UCST, in which the Wilson loop could be evaluated exactly, was 
discussed in Ref. [1] and investigated in the 4D case in Ref. [4], where the string tension of the 
Nambu-Goto term and the coupling constant of the rigidity term were calculated. The latter 
one occured to be negative, which means that the obtained (Euclidean) string effective action is 
stable [5]. Also the 4D UCST action was derived in Ref. [4] by performing the exact duality 
transformation, while in Ref. [1] it was done only at the semiclassical level. 

However, as it was first pointed out already in Ref. [6], the non- Abelian generalization of the 
above described antisymmetric tensor theory is quite difficult. 

Another, more phenomenological, but on the other hand adapted to the non-Abelian case 
approach to the problem of the string representation of the confining phase of gauge theories 
was developed in Refs. [7-9]. It is based on the Method of Vacuum Correlators [10] (MVC) 
and allows one to get the information about the gluodynamics string effective action from the 
expansion of the Wilson average, considered as a statistical weight of this string theory, in powers 
of the correlation length of the vacuum. In this way, in Ref. [7] the sign of the rigidity term 
was found to be negative according to the present lattice data, which means that the MVC is in 
agreement with the dual superconductor model of confinement due to Ref. [5]. In Ref. [8], the 
string effective action obtained in Ref. [7], was applied to the derivation of the correction to the 
Hamiltonian of the QCD string with quarks [11] due to the rigidity term. However, up to now 
it is not clear how the summation over the world-sheets could appear within this approach. The 
most appropriate way in this direction might lie in the accounting for the perturbative gluons, 
which propagate inside the Wilson loop and, as it was argued in Ref. [12], should generate world- 
sheet excitations. In Ref. [9], this conjecture was elaborated out by virtue of the integration 
over perturbative gluons in the Wilson average, which in the lowest order of perturbation theory 
leads to the interaction of the world-sheet elements via the exchanges of perturbative gluons, 
propagating in the nonperturbative background. Expanding this interaction in powers of the 
derivatives w.r.t. the world-sheet coordinates, we finally get in the lowest order of this so-called 
curvature expansion some definite correction to the rigid string coupling constant, while the string 
tension of the Nambu-Goto term acquires no additional contribution due to perturbative gluons 
and keeps its pure nonperturbative value. 
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In the present Letter, we shall propose an alternative method of description of fluctuations 
of the gluodynamics string world-sheet. To this end, we shall first demonstrate in the next 
Section that the gluodynamics string emerging within the bilocal approximation to the MVC 
could be effectively described by the same action of the massive Kalb-Ramond field interacting 
with the string, which describes the low-energy limit of the UCST. This observation will be 
proved by the calculation of the coefficient functions standing at two Kronecker structures (which 
describe surface and boundary terms in the string effective action) in the propagator of the Kalb- 
Ramond field, and further comparison of them with the values of the coefficient functions standing 
at the corresponding structures in the gauge-invariant correlator of two gluonic field strength 
tensors. The latter ones are unfortunately yet not found from the gluodynamics Lagrangian, 
while equations for correlators were obtained in Refs. [13] and [14] and then investigated in Ref. 
[15] by making use of the stochastic quantization method (see also Ref. [16], where alternative 
equations for correlators following from the non-Abelian Bianchi identities, which were proposed 
in Ref. [17] and generalized in Ref. [13], were investigated). That is why, we shall compare the 
coefficient functions, standing in the propagator of the Kalb-Ramond field, with the lattice data 
[18] concerning the behaviour of the corresponding coefficient functions in gluodynamics (see also 
Ref. [19], where these functions were measured in QCD with dynamical fermions). We shall 
see that the mass of the Kalb-Ramond field plays the role of the inverse correlation length of 
the vacuum, so that in the sum rules' limit, when the correlation length of the vacuum tends to 
infinity (see the last Ref. in [10]), which in our model corresponds to the case of the massless 
Kalb-Ramond field, the string picture disappears, and we are left with the boundary terms only. 
The points described above will be the topics of the next Section. 

In Section 3, we shall use the model of the gluodynamics string proposed in Section 2, to 
describe string world-sheet excitations. Correspondingly, for the low-energy limit of the UCST this 
will be an exact procedure rather than a model dependent approach. The world-sheet excitations 
will be described with the help of the background field method developed in Ref. [20] for nonlinear 
sigma models. Namely, we shall split the world-sheet coordinate into a background and quantum 
fluctuations, after which upon the integration over the Kalb-Ramond field we shall derive for the 
latter ones a quadratic action, which contains several nontrivial couplings of quantum fluctuations 
with the background world-sheet with- and without derivatives. 

The main results of the Letter are summarized in the Conclusion. 

In the Appendix, we perform rather nontrivial integration over the Kalb-Ramond field in the 
expression for the 4D UCST partition function. 

2. A Unified Description of the Gluodynamics String and the Low- 
Energy Limit of the UCST 

The partition function of the 4D UCST (which is nothing else but the Wilson average in the 
Euclidean 4D compact QED) in the low-energy limit has the form [1,4] 

(W(C)) = N J DB^exp 

Here 

is a strength tensor of the field B^, 



I 
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T^(x) = J d<jp U (x(£))5(x - x(£)) 



is the vorticity tensor current, e is a dimensionless coupling constant, and A = ^f-^/z, where 

const . 

z ~ e ^ , and Ao is a cutoff which is necessary in 4D. It is worth mentioning, that historically 
Eq. (1) has been for the first time considered in Refs. [21] and [22]. In Ref. [21], such a partition 
function has been introduced for the purposes of construction of the confining medium based on 
an antisymmetric tensor order parameter, whereas in Ref. [22] it appeared during formulation of 
a local quantum field theory of charges and monopoles. 

Gaussian integration in (1) is carried out in the Appendix and yields the UCST low-energy 
action in the form 



Here 



where 



Sucst = J da Xv {x) j da w (x') (B X u(x)Bpp(x')} . 
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In Eqs. (4) and (5), m = j is the mass of the Kalb-Ramond field, i-Q's, i — 0, 1,2,3, stand for 
the Macdonald functions, and one can show that the term 

J da Xv {x) J da w {x') (Bx^B^x'))^ 

could be rewritten as a boundary one, due to which one can immediately establish a correspondence 
between (W(C)) = exp (— SVjcst) an d the Wilson average written within the bilocal approximation 
to the MVC This correspondence yields the following values of the coefficient functions D and 
Di, which parametrize the gauge-invariant correlator of two gluonic field strength tensors 
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e 2 m 
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K 1 (m |x|) Tn ( . ... N 



(7) 



Notice, that the derivative (or curvature) expansion of the action (2), which in this case is 
equivalent to the ^-expansion, was performed in Ref. [4] and up to a common constant positive 
factor has the form 



S UCST = m 2 K J d 2 CV9- \ J d 2 CV99 ab (^V) (^V) + (^) / ds\ 



dx^ dx^ . . 

^7^7- (8) 



Here and ^ stand for the induced metric and the extrinsic curvature tensor of the string 
world-sheet respectively, 



+ °°dt 



f(y) = J T Mt), 



y 

x M (s) in the last term on the R.H.S. of Eq. (8) parametrizes the contour C in Eq. (1), x M (0) = 
x M (l), and we have omitted the full derivative term of the form / d 2 ^^/gR, where R is a scalar 
curvature of the world-sheet. As it was already mentioned in the Introduction, the coupling 
constant of the rigidity term in this expansion is negative, which means the stability of the string. 

Asymptotic behaviours of the functions D and D 1 following from Eqs. (6) and (7) at \x\ <C — 
and \x\ ^> — read 

^ e 2 m 2 . , 

D ~ > 7T5' 9 

Di — > < ^—s (10) 

2tt 2 (x 2 ) 2 V ; 

and 

8V27T2 ( m | x |)2 

2 m 4 -m\x\ 

Di — > ■= 3 5- (12) 

4V27T2 (m|x|) 2 

respectively. 

Let us comment on these asymptotic behaviours. Firstly, one can see that Eq. (10) is in 
agreement with the MVC, where due to Ref. [10] at the distances which are much smaller than 
the correlation length of the vacuum T g , 



Di — -P^f- (13) 



16a s (x 2 ) 
3tt {x 2 f 

The difference in the numerical constants in the asymptotic behaviours (10) and (13) is due to the 
colour factor in the one-gluon exchange diagram, which contributes into Eq. (13). This factor, 
which is absent in the asymptotics (10) of the Abelian propagator, could be accounted for by the 
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proper tuning of the charge e of the Kalb-Ramond field, if we replace the running a s (x 2 ) in Eq. 
(13) by some fixed value, say approximate a s (x 2 ) by its "frozen" value, which it acquires at the 
confinement scale [23] where a is the string tension of the Nambu-Goto term. 

However Eq. (9) agrees with MVC only in the lowest order of perturbation theory, when the 
perturbative part of the function D is absent, and D tends to the gluonic condensate at \x\ — > 0, 
i.e. this agreement is only in a sense that DcDi. One can see that Eq. (9) does not coincide 
with the behaviour of the function D (x 2 ) of the type 



(* 2 ) 2 



In (mV) + 



where a, (3 stand for some constants, and M is a certain mass parameter, which takes place in the 
next-to- leading order of perturbation theory [24]. This is due to the fact that such a behaviour is 
a specific property of non-Abelian theories [25] and presumably could not be reproduced by any 
Abelian model. 

Secondly, comparing asymptotic behaviours (11) and (12) we see that the function D 1 falls off 
much faster than the function D, which is in agreement with the lattice data [18]. The exponential 
falls-off of the functions D and D\ also agree with Ref. [18], while the preexponential power-like 
behaviours do not. However the T^-fits of these behaviours used in Refs. [18] and [19] were 

motivated by the short-distance asymptotics of the function D ± , which according to Eq. (10) is 
reproduced in our model as well. 

Eqs. (11) and (12) tell us once more that the mass m of the Kalb-Ramond field in our approach 
should be treated as an inverse correlation length of the vacuum T g , so that in the "string limit 
of QCD" [26], when the string tension of the Nambu-Goto term, a, is kept fixed at vanishing 



T g , we have a correspondence of the type m ~ y I n the opposite regime when the Kalb- 
Ramond field is massless (or equivalently, in the strong coupling regime of the theory (1)), which 
corresponds to the QCD sum rules' case (see the last Ref. in [10]), the "confining" part (4) of the 
propagator of the Kalb-Ramond field vanishes, and we are left only with the "boundary" part (5), 
whose contribution to Sucst takes the form 



(x — x') 



.A 2 ' 



J da Xv (x) j da w (x') (Bx^B^x'))^ = l^f dx »f dx '^\ 

c c 

which could be anticipated from the very beginning. This result is also in agreement with Ref. 
[4], where it was shown that in the strong coupling regime of the theory (1) strings are completely 
suppressed when A — > +oo. 

3. Description of the World-Sheet Excitations 

In order to describe fluctuations of strings in the model (1), we shall adopt the background- 
field method developed in Ref. [20] for the nonlinear sigma models. The three differences of our 
case with the one considered in Ref. [20] are the absence of the Polyakov term in the action, 
flatness of the field manifold (i.e. the space-time), since we consider an arbitrary Wilson loop, not 
necessarily lieing on the unit sphere, and the necessity of the eventual integration over the field 
B^. Consequently, the geodesies passing through the background manifold y M (£) and the excited 
manifold x M (£) = y + £/x(£)> where stands for the world-sheet fluctuation, are straight, 

Pfi(£,s) = y M (0 + sz (t(Q, where s denotes the arc-length parameter, < s < 1. The expansion 
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of the action (2) in powers of quantum fluctuations z^(0 could be performed by virtue of the 
following generating functional, which is the arc-dependent term describing the interaction of the 
Kalb-Ramond field with the string in Eq. (1), 



I [p(£, s)]=if d\B, u [p(£, s)\ e ab (d aP ^ s)) (d bPu (Z, s)) . 
Then the term containing n quantum fluctuations reads as 

1 d n 



j(n) 



n\ ds' 1 



s=0 



and we get from Eq. (14) 



= i J da, u (y(0)B,u [y(0] , 

/« = i J daMO>x(OH„x [y(0] , 



and 



(14) 



(15) 
(16) 



J( 2 ) = i J d^z u (Oe ab (d a y^)) H u , x [y(0\ + ^(0 (d b y x (0) d a H vllX [y(0\j , (17) 

where during the derivation of Eqs. (16) and (17) we have omitted several full derivative terms. 

It is worth mentioning, that as it was discussed in Ref. [20], the terms (15)-(17) are necessary 
and sufficient to determine all one-loop ultraviolet divergences in the theory (14) at s = 0. This 
statement is important for the model of the gluodynamics string with the action (2), proposed in 
the previous Section, since in the UCST case, Eq. (2) describes only its low-energy effective action, 
which does not suffer from the ultraviolet divergences. It is also important for the Abelian Higgs 
Model in the London limit [27], where Eq. (1) with an additional integration over world-sheets is 
nothing else but the expression for the 't Hooft loop average. 

In order to get the desirable action, quadratic in quantum fluctuations z^), we shall first 
carry out the integral 

N J DB, V exp \- J dx y^ 2 Hl„ x + ±B%)j + /<°> + 
It occurs to be equal to 



exp 



doxv (y(0) I da, p (y(0) ( (B Xv [y(0] B„ W)]f ] + 2^(0 



d 



dy a (0 



(B,AymB» P [y(a) {1) +z«(OMe 



d 2 



'M\W 



T AB Xu W)\B, P W)\) 



(18) 



where from now on we shall omit all the boundary terms. 

Secondly, one should substitute the saddle-point of the integral 
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into Eq. (17). This saddle-point reads 



12 .\2 H Z>* + Ae 2 B l>) +/(0) 



■ i 3 
xe m 



2jr2 j \y(()-y((')\ • 

and upon its substitution into Eq. (17), accounting for Eq. (18), and making use of Eq. (4), we 
finally get the following value of the action quadratic in quantum fluctuations 



quadr. 



e 2 m 3 
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(19) 



where 



/ ^ (y(O) (y(0 - v(0)„ + / *v (v(O) (y(0 - v(0)a + / ^ (y(O) (y(0 - v(0)„ 



— ^-- + ^(Ko + ^) 

|y(0-y(C)l 2 



(20) 



(y(0-y(£')) 2 

and everywhere in Eqs. (19) and (20) the arguments of the Macdonald functions are the same, 
m \y(0 ~ y(£')l- We have also kept the pure background term on the R.H.S. of Eq. (19). The 
other terms yield the couplings of quantum fluctuations with the background world-sheet. All of 
these terms except one do not contain the derivatives of quantum fluctuations. 



4. Conclusion 

In the present Letter, we have modelled the gluodynamics string effective action by the theory 
of the massive Abelian Kalb-Ramond field, interacting with the string. The partition function 
(1) of this theory is nothing else but the low-energy expression for the partition function of the 
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UCST, which is in fact the low-energy limit of the Wilson average in the 4D compact QED, 
rewritten in terms of the dual antisymmetric tensor field. This observation provides us with a 
unified description of the gluodynamics string and the UCST. 

An approach to the gluodynamics string with the help of the theory (1) has been justified 
in Section 2, where it has been shown that the small- and large-distance asymptotic behaviours 
of the coefficient functions, which parametrize the gauge-invariant correlator of two gluonic field 
strength tensors within the MVC and could be extracted from the lattice measurements, are in a 
good agreement with the ones of the corresponding functions standing in the propagator of the 
Kalb-Ramond field. It has also been demonstrated that the mass of the Kalb-Ramond field plays 
the role of the inverse correlation length of the vacuum, so that in the massless limit strings are 
suppressed, which is in agreement with Refs. [4] and [26]. 

In Section 3, we have applied the background field formalism, developed in Ref. [20] for the 
nonlinear sigma models, in order to describe fluctuations of strings in our model. To this end, 
we have splitted the string world-sheet coordinate into a background part, corresponding to some 
given world-sheet (e.g. one with the minimal area), and quantum fluctuations, and using the arc- 
dependent string action as a generating functional, performed the expansion of the term, which 
describes the interaction of the string with the Kalb-Ramond field, up to the second order in 
quantum fluctuations. Finally, upon the integration over the Kalb-Ramond field, we have derived 
an action, which is quadratic in quantum fluctuations and contains the pure background part and 
the terms describing the interaction of the background world-sheet with the quantum fluctuations. 
This action is given by formulae (19) and (20) and describes fluctuations of strings in the UCST 
and in our model of the gluodynamics string. 
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Appendix. Integration over the Kalb-Ramond Field in the UCST 
Partition Function (1) 

According to the general rule, in order to calculate the Gaussian integral (1) one should substi- 
tute its saddle-point value back into the integrand. The saddle-point equation in the momentum 
representation reads 

This equation can be most easily solved by rewriting it in the following way 

(p 2 Px„, a p + m 2 lx„, a p) B^ {p) = 2tA 2 T Xu (p), (A.l) 

where we have introduced the following projection operators (see e.g. Appendix to the last Ref. 
in [15]; functional generalization of these operators has been introduced in Ref. [3]) 
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and 



where V u 



1 



8n.v — ^r 1 - These projection operators possess the following properties 

l/i^,Ap ^-up,,Xp lp^,pA ~^-\p,p,V7 



l/i^,AplAp,a/3 ^-p,i/,af3 



(the same properties hold for P^^), and 

P^,Ap (1 - P, 



\p,aj3 



0. 



By virtue of properties (A.2)-(A.4), the solution of Eq. (A.l) reads 



2tA 2 



pi _|_ m 2 



1 + ^(1-P) 

77T 



\v,af3 



(A2) 



(A3) 



(A4) 



which, once being substituted back into partition function (1), yields for it the following expression 



(W(C)) = exp 



-A 



: f dp 1 
J (2tt)V + 



m z 



P 



1 + S(1-P) 



rrr 



T^{-p)T a p(p) 



- p,v,af3 



(A5) 



Rewriting Eq. (A. 5) in the coordinate representation we arrive at Eq. (2). From Eq. (A. 5) one 
can easily see that the term on its R.H.S. proportional to the projection operator 1 — P yields in 
the coordinate representation the boundary term, i.e. Eq. (5), whereas the term proportional to 
the projection operator 1 yields Eq. (4). 
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